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ABSTRACT 



The hypothesis that the cyanomethyl anion CIUCN is responsible for the relatively narrow diffuse interstellar band (DIB) at 8037.8 + 
0. 1 5 A is examined with reference to new observational data. The 0[j absorption band arising from the 1 B i - X 1 A' transition from the 
electronic ground state to the first dipole-bound state of the anion is calculated for a rotational temperature of 2.7 K using literature 
spectroscopic parameters and results in a rotational contour with a peak wavelength of 8037.78 A. By comparison with diffuse band 
and atomic line absorption spectra of eight heavily-reddened Galactic sightlines, CIUCN" is found to be a plausible carrier of the 
/18037 diffuse interstellar band provided the rotational contour is Doppler-broadened with a b parameter between 16 and 33 km s 
that depends on the specific sightline. Convolution of the calculated ClUCN" transitions with the optical depth profile of interstellar 
Tin results in a good match with the profile of the narrow /18037 DIB observed towards HD 183143, HD 168112 and Cyg OB2 8a. 
The rotational level populations may be influenced by nuclear spin statistics, resulting in the appearance of additional transitions from 
K a = 1 of ortho CHiCN~ near 8025 and 8050 A that are not seen in currently available interstellar spectra. For CH^CN" to be the 
carrier of the /18037 diffuse interstellar band, either a) there must be mechanisms that convert CUCN" from the ortho to the para 
form, or b) the chemistry that forms CH2CN~ must result in a population of K' a ' levels approaching a Boltzmann distribution near 3 K. 

Key words. Astrochemistry — ISM: lines and bands - ISM: molecules - ISM: atoms - ISM: clouds 



1. Introduction 

The origin of the unidentified diffuse interstellar bands (DIBs) 
remains one of the greatest challenges in astronomical s pec- 
troscopy. Th e subject has been reviewed by iHerbigl d 1995b and 
ISarrd d2006l) who have highlighted research that points towards 
organic molecules as likely candid ates for at least some of the 
more than 300 DIBs. ISarrd d2000l) has put forward a hypothe- 
sis that 'some, possibly many' of the diffuse interstellar bands 
arise from electronic transitions between ground and dipole- 
bound states of negatively charged polar molecules or small po- 
lar grains. It was noted that the r Qo(l) line of the origin band of 
the l Bi - X 'A' transition of CH 2 CN occurs at 8037.8 A, in 
good correspondence with the peak absorption wavelength of a 
diffus e interstellar band at 8037.9 + 0.3 A dGalazutdinov et al.L 
2000). 



2. Molecular anions as diffuse band carriers 

2.1. Previous studies 

The electronic absorption s p ectrum of in the visible region 
was studied by iTulei et at] dl998l) and five optical absorption 
bands were reported to match with the wa velengths of known 
DIBs. Models of diffuse cloud chemistry by iRufHe etaTI d 19991) 
that incorporated the desorption of seed molecules from grain 
surfaces were able to reproduce large abundances of (approx- 
imately 10~ 9 «hX an d also of C7H , albeit only under a limited 



set of physical and chemical conditions and only for short peri- 
ods of time. The high electron affinity and density of vibrational 
states of C7 permits rapid radiative election attachment such 
that anion formation has b een calculated to occur ev ery time an 
electron collides with C 7 dTerzieva & Herbstl |2000|) . However, 
the early promise of as a DIB carrier was quashed when 
high resolution ob s ervati o nal and laborato r y spe c troscopy by 
iGalazutdinov et al.l dl999]) ISarre & Kendalll (12000b . lLakin et al.l 
d2000b and lMcCalletalTd2001h identified that the match be- 
tween the DIBs and the wavelengths, strengths and profiles of 
the optical absorption bands was too poor to constitute an 
assignment. 

Based on labo ratory spectroscopy it was suggested by 
iGuthe et all d2001b that CH 2 CC might be a diffuse b and car- 
rier but a detailed assessment by iMcCall et al.l d2002bl) showed 
that the wavelength match to the diffuse band at 6993 A was 
not acceptable and commented that the non-observation of ad- 
ditional K sub-bands argued against an assignment to CH2CC~. 
The CH2CN molecule has many spectroscopic characteristics 
in common with CH2CC , both molecules being near-prolate 
asymmetric tops with ortho and para forms and both possessing 
dipole-bound excited states, but CH2CN has a closed electronic 
shell. 



2.2. Anion chemistry in the ISM 



(with kinetic temperatures of 
3 and strong UV radiation 
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Chemistry in the diffuse ISM 
~ 50 K, densities ~ 100 cm~ D and strong 
fields) has been shown to be rich and complex by observations 
whi ch include the detec tion o f polyatomic mole cules such as 
H+ dMcCall et all 11998b . C 3 dMaier et all 1200 ll) . and HCO + , 
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HCN, H 2 CO C 2 H an d c-C 3 H 2 dLucas & LiszlL fl996allbl l2000t 
iLiszt & Lucasl 1200 lb at high fractional abundances similar to 
those observed in dark clouds. Conventional models of gas- 
phase diffuse cloud chemistry struggle to explain how such 
high molecular abundances can be maintaine d in the interstellar 
UV ra diati on field. However, as disc ussed bv lDulev & Williamsl 
dl984l) and lHall & Williamsl ( 1 1995b . the availability of erosion 
products from carbonaceous dust grains, including polyynes (C„ 
chains) and small molecules such as diacetylene (C4H2), could 
provide a reservoir of reagents to fuel a complex network of or- 
ganic reactions in diffuse clouds and other photon-dominated re- 
gions. Large molecules should be able to resist photodissociation 
in the interstellar UV field by undergoing internal conversion in 
which the energy of absorbed photons is redistributed into the 
vibrational modes of the molecule due to the high density of 
states. To date the role of anions in interstellar reaction networks 
has generally been considered to be minimal although some con- 
sideration of the ir possible d etection in the interstellar medium 
has been made ( ISarreL Il980h . No reaction pathways involving 
anions were present in the 'compreh ensive' diffuse cloud mod- 
els of Ivan Dishoeck & Black! dl986b (where ~ 500 reactions 
were modelled), and only a few atomic and diatomic anions are 
present in the UMIST 1999 database of ~ 4000 astrochemicalre - 
actions involving ~ 400 different species dLe Teuff et al. 
This is at least partly due to the lack of observational evidence 
for anions in space, but also reflects a common view that anions 
should be rapidly photoionised in the interstellar UV field. 

The abundance of interstellar anions in relatively diffuse me- 
dia depends in la rge part on the balance betw een radiative elec- 
tron attachment dDalgarno & McCravl 1 1973b and photo detach- 
ment b y th e interstellar rad i ation field. As discussed by iHerbslI 
d!98lb and lHerbst & Petrid dl997b . a viable route for molecular 
anion formation in the ISM involves attachment of a free elec- 
tron via the formation of an excited temporary anion that sta- 
bilises through radiative transitions to the ground state. Provided 
the radiative transition rate is high, anions should be able to form 
at a sufficient rate that - relative to their neutral counterparts - 
appreciable anion abundances may occur depending on the local 
interstellar gas pressure and radiation field. Chemical models of 
the crrcums t ellar e nvelope of the carbon star IRC+10°216 by 
iMillar et"aT1 d2000b predict observable abundances of C„~, C„H 
and C„N~. The discovery of the microwave signatur e of C(,H in 
IRC+10°216 and TMC-1 bv lMcCarthv et alj [2006b at fractional 
abundances (relative to the neutral species CeH) of about 2.5% 
and 1 %, respectively, proves that anions arise in some regions of 
space in significant quantities. 

2.3. Dipole-bound states 

Anions derived from strongly polar neutral parent molecules 
(or grains) can support one or more 'dipole-bound' electronic 
states at energies near to the detachment threshold and can po- 
tentially form the excited states in electronic absorption spec- 
tra in the visible region of the spectru m. Dipole-bound elec- 
tronic states were first predicted by iFermi & Teller! d 1947b 
who showed that a theoretical point dipole with a sufficiently 
large dipole moment (> 1.625 D) can bind an elec t ron in a 
diffuse orbital. Theoretical studies dCrawfordl Il970t iGarrettl 
ll978HGutsev & Adamowiczlll995allbl|ch and experimental work 
dMoran et all Il987t iBrinkman et all Il995l iLvkke etall. Il987t 
iDesfrancois et all 1 1996b has since demonstrated that in practice, 
dipole moments fi > 2.5 D are required for a molecule to sup- 
port at least one dipo le-bound electronic st ate. Laboratory stud- 
ies of nitromethane dCompton et all 1 1996b and of cyanoacety- 



Target 


Sp. type 


V 


Eb-v 


d (pc) 


Cyg OB2 5 


07 Ie 


9.2 


1.94 


700 


Cyg OB2 8a 


O6 I 


9.0 


1.59 


1000 


Cyg OB2 12 


n ^ 1 n t 

B2I-B8Ia 


1 1.4 


2.8-3.4 


600 


T Tl~"\ 1 f O 1 1 O 

HD 168112 


05 


8.6 


1.01 


2000 


1 1 \J 1 l)7Ujt 


It 

DJ la 


o._ 


1 -1 
1 


1 1 on 

1 1UU 


HD 183143 


B7 la 


6.9 


1.24 


650 


HD 186745 


B8Ia 


7.1 


0.96 


1600 


HD 229196 


05 


8.6 


1.22 


1100 


iCyg 


B5 V 


4.6 


0.04 


50 


f Peg 


B8.5 V 


3.4 


0.00 


50 



Table 1. Table showing sightline data for /18037 observa- 
tions. For the Eb-v v alues, B - V photometry are from 
iPerryman & ESAl d 19971) and spectral types from the modal av- 
erages of all data referenced in the SIMBAD database (URL: 
http://simbad.u-strasbg.fr/sim-fid.pl); intrinsic stellar photome- 
try is from Wegner ( 1994) with the except ion of Cyg OB 2 5 and 
8a for which photometry was taken from iMassey & ThompsonI 
dl99lb and Cyg OB2 12 which is of uncertain s pectra l 
type and photometry as discussed by iGredel & Munchl d 19941) . 
Approximate stellar distances d (accurate to around +50%), 
w ere calcu l ated fr om the spectral types and absolute magnitudes 
of lWegnerl d2006l) (stellar apparent magnitudes were dereddened 
assuming a ratio of visual to selective extinction (Ry) of 3.1). 



lene and uracil dSommerfeldl 12005b indicate that the presence 
of a near-threshold dipole-bound state can act as a 'doorway' 
for electron capture, and that strong electronic transitions occur 
from the dipole-bound state to the ground state, resulting in rapid 
radiative stabilisation of the anion. It is possible that this effect 
may enhance the electron capture rate for the strongly dipolar 
CeH molecule and t herefore contribute tow ards the abundance 
of C 6 H observed bv lMcCarthv et al.1 (12006b . 

3. Spectroscopic observations and data reduction 

Medium resolution optical spectroscopic observations of early- 
type Galactic stars were performed using the High Resolution 
Echelle Spectrograph (HIRES) of the W. M. Keck Observatory, 
Hawaii by G. H. Herbig between 1995 and 1997 (private com- 
munication), who kindly made available the raw science expo- 
sures of eight heavily-reddened Galactic sightlines (shown in 
Table[T|i, including flat fields, Th/Ar and bias frames. Reductions 
were carried out using standard iraf echelle routines. In the re- 
gion of the /18037 DIB special care was taken to eliminate tel- 
luric lines and fringing residuals by division with high S/N stan- 
dard star spectra. The S/N of the reduced spectra is typically be- 
tween 400 and 800 per pixel, with a resolving power of 42 500 
(measured from the average of several unblended Th/Ar lines 
over the wavelength range of ~ 7000 to 10000 A). In addition, 
near-UV spectra from ~ 3200 to 4000 A were obtained for the 
sight-lines towards Cyg OB2 8a, HD 168112 and HD 183143 
and include spectra of a range of important interstellar species 
including Ca 11, Ti 11 and various diatomic molecules that provide 
valuable information on the velocity distribution of the interstel- 
lar gases. The near-UV spectra typically have S/N ~ 300 with a 
resolving power of 52 500. Around 8040 A the wavelength drift 
between successive arcs was 0.02 A, and for the near-UV ob- 
servations (around 3600 A) the drift was 0.005 A. Thus, abso- 
lute wavelength calibration errors should be AA/A < 2.5 x 10" 6 
(0.75 km s- 1 ) for the red/NIR and < 1.4 x 10~ 6 (0.4 km s" 1 ) for 
the near-UV spectra. 
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4. Spectroscopy of CH 2 CN 



The neutral cyanomethyl radical CH2CN has a dip ole mo- 
ment between c. 3.5 and 4.0 D dOzeki et al l 12004 and an 
electron affinity of c. 1.55 eV dLvkke et all Il987l) . The first 
dipole-bound state of the anion l i es ^ 10 meV bel o w the 
tonisation continuum dLvkke et all fl987t iMoran et all fl987t 
iGutsev & Adamowicj.il 995 al) . The CH^CN radical is planar but 
on attachment of an electron the hydrogen atoms move out of 
pla ne, resulting in a X 'A' electron ic ground state in C s symme- 
try dGutsev & Adamowiczl Il995al) with the mirror plane along 
the C-C-N backbone. In the excited 'Bi state, the dipole-bound 
electron is very diffuse with only a weakly perturbative inter- 
action with the rest of the molecule. Thus in the dipole-bound 
state the nuclear framework assumes a geometry which is al- 
most identical to that of the ground-state neutral radical. The 
'Bi - X 'A' transition is of perpendicular type dLvkke et all 
1 1987b with the transition dipole orientated along the c-axis and 
selection rules AJ = +1, AK a = +1 (K a is used although it 
is not strictly a good quantum number). The rotational con- 
stants of CH2CN were determined to very high accuracy by 
iLvkke et all ( fl98l using high-resolution fast-ion-beam autode- 
tachment spectroscopy. The rotational energy level structure and 
som e examples of 'Bi — X 1 A' transitions are shown in Figure 
4. of lLvkke etal] (119871) . 

Nuclear spin statistics dictate that there exist ortho (o) 
and para (p) forms of CH2CN as occurs in the isoelec- 
tronic molecule cyanamide (NH2CN), discussed bv lMillen et alj 
dl962l) . If the ortho and para CHbCN abundances reflect the 
o : p (3 : 1) nuclear spin degeneracies, transitions of ortho 
CH2CN would be three times stronger than equivalent transi- 
tions of the para form. It not clear whether thermalisation of 
the rotational energy levels of CH>CN would take place, but 
based on observations of other interstellar molecules, it is evi- 
dent that ortho-to-para abundance ratios do not necessarily fol- 
low the nuclea r spin degeneracies. For example, according to 
iFlower & Watt! d 1984b . interconversion between ortho and para 
H2 occurs via collisional proton exchange in the reaction 0-H2 
+ H + — > P-H2 + H + + hv, which allows the population of the 
lowest J levels to reach thermal equilibrium with the ISM such 
that o : p is skewed away from the value of 3 : 1 . In the case of 
CH2CN the photodetachment (destruction) rate has been calcu- 
lated to be approximately 1 .2x 10~ 8 s -1 (E. Herbst & T. J. Millar, 
private communication); given the relatively slow rate of conver- 
sion between gas-phase orth o and para H2 by prot on exchange 



5. Results 



with H + (3 x 10 10 n H + s" 1 : IFlower & Wattl J 1984 ). it follows 



that in the neutral ISM where H + densities are generally less 
than 0.1 cm 3 , collisional exchange with gas-phase CF^CN - is 
unlikely to alter significantly the ortho-to-para population ratio. 

For molecules with equivalent H atoms the rotational level 
populations (and o : p ratios) can give clues as to the chemical 
reactions in which the molecule participates. For example, if a 
molecule is formed by reaction with H2, then the o : p ratio of 
the product may be influenced by the o : p ratio of the reage nt 
H2 (see the case of 0C3H2, studied bylTakakuwa et all d2001l) ). 
Also, according to iDickens & Irvine! d 1999b . the observed o : 
p ratio of H2CO, (which has the same nuclear spin statistical 
properties as CH2CN ), indicates that the molecules probably 
formed in thermal equilibrium with cold dust grains. Both of 
these types of reactions could potentially influence the o : p 
ratio ofCH 2 CN . 



Using the molecular constants of Lvkke et alj dl987b and the 
ASYROT FORTRAN code dBirss & RamsavL 11984 . the % - X 'A' 
transitions of CH2CN were computed for a distribution of ro- 
tational level populations in equilibrium with the 2. 74 K CMB. 
Isoele ctronic with NH2CN (dipole moment 4.32 D; iTvler et all 
dl972f) ). the cyanomethyl anion is expected to have a large dipole 
moment (~ 4 D) such that in the diffuse ISM the radiative tran- 
sition rate is sufficient that collisional excitation of the molecule 
should be small. Small changes in the level of (J) rotational ex- 
citation of the molecule in fact have relatively little impact on 
the spectrum, which is dominated by the Q branch of the tran- 
sition at 8037.8 A. Convolved with a Gaussian with Doppler 
b — 1 km s _1 and at a resolving power of 42 500, the calcu- 
lated spectrum is displayed in Figure[T]for comparison with the 
observed HIRES spectra. 

Examination of the spectrum of the standard star A Cyg 
shows no evidence for significant stellar lines, fringing arti- 
facts or telluric residuals across the wavelength region. The 
most prominent features are the /18026 DIB, which is among 
the narrowest known diffuse bands, and an absorption feature 
around 8037 A. The ,18037 DIB is relatively narrow (FWHM ~ 
1.3-2 A), but is overlapped by a broader absorption centered 
at around 8040 A (noted bv lHerbig & Lekal 1 19911) that is most 
prominent towards Cyg OB2 12 and HD 183143. The ,18040 
component has a FWHM of around 4 A and forms an absorption 
peak at 8040.7 + 0.3 A. In the co-added (Ki rest frame) spectra 
(shown in Figure|2), Gaussian fits to the peak of the /18037 DIB 
give a wavelength of 8037.8 + 0.15 A for this diffuse interstellar 
band. 

The peak wavelength match of /18037 with the 2.7 K 'Bi - 
X 'A' CH2CN spectrum shown in Figure Q] is very good. 
However, the computed spectrum clearly contains significant 
fine structure that is not present in any of the observed ,18037 
profiles. The Q branch of the K'J = transition creates the 
prominent peak at 8037.78 A. The P branch is of rather low 
intensity compared to the R branch that produces the set of 
lines between 8035 and 8037.2 A. Evidence for asymmetry of 
,18037 can be seen in the blue degradation of the DIB spectra 
for the three Cyg OB2 sightlines and HD 183143, HD 186745 
and HD 229196. The band shows significant evidence for profile 
variability across all of the sightlines, especially in comparing 
HD 186745 and HD 183143 with Cyg OB2. 

The variable strength and profile of /18040 with respect to 
,18037 suggests that it is probably caused by a different car- 
rier. Towards Cyg OB2 5 and HD 168112, the ,18040 DIB has 
FWHM ~ 1 1 .5 A - much broader than in the other sightlines 
- and perhaps indicates the presence of another broad DIB. To 
isolate the ,18037 DIB for analysis is non-trivial among these 
other contaminating features. However, it is relatively narrow 
which assists its rejection in continuum fitting algorithms such 
that consistently repeatable continuum fits were possible. 

5.1. Ortho and para CH 2 CN- 

As detailed in SectionH] CH 2 CN is expected to exist in the ISM 
in para and ortho forms, with the lowest occupied rotational lev- 
els being J" = 0, K'J = and /" = 1, K'J = 1, respectively. 
The K„ - 1 < — transitions peak around 8037.8 A, whereas 
the K a — < — 1 and the K a — 2 < — 1 transitions peak at 
8049.6 and 8024.8 A respectively. The ,18037 co-added spec- 
trum for all eight sightlines is plotted in Figure [2] with para + 
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l ' 1 ' r 

CH 2 CN (2.7 K) 




0.64 • i -^fiA/t Vr jYw^ 



8025 8030 8035 8040 8045 8050 
K I rest wavelength (A) 

Fig. 1. Normalised (second order Chebyshev polynomial), 
telluric-corrected HIRES spectra of the region around 8037 A. 
Spectra have been Doppler shifted to place the mean K i wave- 
length at rest (using ,l K i rest = 7698.9645 A ( iMortonl [2003h . A 
Cyg is shown as an unreddened standard. Fitted continua are 
shown as dashed lines. The 2.74 K CH 2 CN % - X 'A' ori- 
gin band absorption spectrum is plotted at the top, calculated as- 
suming a single cloud at rest with Doppler b — 1 km s and 
convolved to the resolving power R - 42 500 of the HIRES 
spectra. The KJ = transitions lie between 8035 and 8040 A 
and at this temperature the KJ = 1 transitions at around 8025 
and 8046 A are almost invisibly weak. The absorption feature at 
around 8026 A is an unrelated diffuse interstellar band. 



ortho CH>CN 'Bi - X 'A' transitions, convolved with a rest 
cloud model (vlsr = km s _1 ) with Doppler b — 18 km s _1 and 
an arbitrary intensity scaling. Two KJ population scenarios are 
plotted: (1) a 3 : 1 ratio of odd : even KJ level populations and 
(2) a 2.74 K Boltzmann distribution. 




1.06 



1.04 



1.02 



0.98 



0.96 - 



0.94 

8020 8025 8030 8035 8040 8045 8050 8055 
K I rest wavelength (A) 

Fig. 2. Comparison of computed and observed (averaged) spec- 
tra in the 8037 A region. The bottom trace shows the co- 
added, normalised, telluric-corrected HIRES spectrum around 
,18037 observed towards the stars listed in Table Q] Spectra were 
Doppler-shifted to place the weighted-mean Ki wavelength at 
rest before co-addition. The % - X 'A', \K a = 2 ^- 1, 
K a — 1 < — 0, K a = < — 1} modelled absorption spec- 
trum is shown, and includes the transitions arising from both 
KJ = (para) and KJ = 1 (ortho) CH 2 CN . The middle spec- 
trum was calculated assuming a thermal Boltzmann distribution 
of K a level populations. The upper spectrum shows the result 
taking the KJ populations to be determined solely by a 3 : 1 
nuclear spin-statistical weight ortho : para ratio for odd : even 
KJ levels. The interstellar CH2CN distribution is modelled as 
a single cloud at rest (v = km s _1 ) with Gaussian line-shape 
(Doppler b = 18 kms 1 ), convolved to the resolving power 
R — 42 500 of the HIRES spectra. The respective peak absorp- 
tion wavelengths of the three Doppler-broadened K a sub-bands 
are 8024.8, 8037.7, and 8049.6 A, shown as dotted lines and cor- 
respond to ortho, para and ortho CH2CN respectively. 



The predicted K a = 2 < — 1 transitions (with peak absorp- 
tion at ~ 8024.8 A) partially overlap the narrow ,18026 DIB. 
There is some evidence for a small blue shoulder on this DIB, at 
a wavelength slightly blue of the peak CH2CN K a — 2 < — 1 
absorption wavelength. The K a = < — 1 feature (with peak ab- 
sorption at ~ 8049.6 A) falls close to a small absorption feature 
in the co-added spectrum. In this spectrum, the central depth of 
the A8037 feature is 1.35 ± 0.16% and the central depth of the 
^8049 feature is 0.3 + 0.16% (2<x error estimates derived from 
the RMS noise of the continuum). The ratio of the central depth 
of the /18037 to ,18049 features is consistent with a o : p ratio of 
~ 1 : 2.3, but given the mismatch between the calculated and ob- 
served profile and peak wavelength, it seems unlikely that ,18049 
is caused by 0-CH2CN . 

From Figure [2] it is clear that if interstellar CH2CN has a 
3 : 1 'statistical' ratio of o : p states, it cannot be the carrier of 
the /18037 DIB due to the lack of the ortho transitions at the ex- 
pected strengths in the observed spectra. The /18026 DIB is too 
red relative to the Q branch (main peak) of the K a — 2 < — 1 
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Fig. 3. Telluric-corrected, normalised DIB ,48037 spectra (thin 
black traces) and least-squares fitted CH 2 CN" 'Bi - X 'A', 
K a = 1 < — models (thick red traces). Models employ a sin- 
gle Gaussian interstellar cloud. Radial velocity shift, Doppler 
width, central depth and an additive continuum offset were free 
parameters in the fits. The DIB rest wavelength was set at Ao = 
8037.78 A for the displayed velocity scale. Best-fitting model 
parameters are shown in Table [2] and the inferred CH 2 CISL in- 
terstellar velocity distributions plotted in Figure [4] (lines labeled 
'Mdl'). 



Sightline 


VLSR 


b 


Ws037 


a c 


°7 


Cyg OB2 5 


13.9 


21.6 


32.5 


0.0020 


0.0014 


Cyg OB2 8a 


12.3 


21.0 


30.2 


0.0024 


0.0016 


Cyg OB2 12 


9.4 


24.0 


38.3 


0.0024 


0.0013 


HD 168112 


29.9 


26.4 


24.5 


0.0018 


0.0021 


HD 169034 


13.6 


33.0 


38.7 


0.0036 


0.0024 


HD 183143 


25.0 


18.0 


37.5 


0.0022 


0.0012 


HD 186745 


19.3 


16.2 


24.3 


0.0014 


0.0013 


HD 229196 


14.1 


21.6 


22.0 


0.0020 


0.0012 



Table 2. Least-squares fit parameters of the /18037 CH 2 CN 
models shown in Figure [3] The mean LSR velocity 
(vlsr / km s _1 ), Doppler cloud width (b / km s _1 ) and equivalent 
width of the model /I8037 DIB (W%o3i I mA) are given along with 
the normalised continuum RMS of the observed spectra (cr c ), 
and the RMS of the model fit residuals (cry). 



transition to constitute a spectroscopic match with the model 
calculation. If however the K a level populations have, or ap- 
proach, a Boltzmann distribution, then the transitions originating 
in K a = 1 are sufficiently weak for CH 2 CN to be a plausible 
carrier of the ,18037 diffuse interstellar band. 




-100 -50 



100 -50 



50 100 



LSR velocity (km s' 



-1, 



Fig. 4. Each of the eight panels above shows LSR velocity- 
space apparent column density profiles for interstellar species for 
which data were available. Data for the K 1 7699 A, CH + 3958 A, 
Can 3968 A and Tin 3384 A lines are shown. Hypothetical 
model CFLGST velocity distributions (labeled 'Mdl') are shown 
at the top of each panel, calculated based on the fits shown in 
Figure [3] (parameters given in Table |2| with the assumption that 
the narrow /18037 DIB is caused by the origin band 'Bi - X 'A' 
transitions of CH 2 CN at 2.74 K, and that the interstellar ve- 
locity distribution of the gas is Gaussian. The interstellar Ki, 
CH + , Can and Tin apparent column densities (N a ) were calcu- 
lated from the normalised HIRES spectra (I/Iq, where Iq is the 
co ntinuum intensity) accordin g to N a ln(/o//) (see for exam- 
ple [Savage]&^e^SacH[i99ll). For the right-hand panel, H i col- 
umn density profiles ar e from the Leiden/ Arg entine/Bonn (LAB) 
Survey of Galactic H i dKalberla et aUl2005l) . 



5.2. Modelling the A8037 diffuse interstellar band 

As shown by Figures Q] and [2] a large Doppler broadening of 
the interstellar gas is required for CH2CN~ to be the carrier of 
/18037. Based on the assumption that the /18037 DIB is caused 
by the origin band % - X 'A', K'J = transitions of CH 2 CN" 
at 2.74 K, it is possible to use the shape of the DIB to in- 
fer the shape of the hypothetical interstellar CH2CN~ velocity 
distribution. This distribution may then be compared with ob- 
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Fig. 5. 2.74 K CH 2 CN % - X 'A', K'J = transitions con- 
volved with the HIRES PSF and the interstellar Tin profiles 
(solid lines), and the Ki profiles (dashed lines), observed to- 
wards HD 183143, HD 168112 and Cyg OB2 8a, overlaid on 
the normalised, telluric-corrected /18037 DIB profiles (thin black 
histograms). 



served interstellar velocity profiles for different species in the 
observed sightlines. The computed 2.74 K absorption spectrum 
was subjected to a least-squares analysis whereby it was con- 
volved with a single Gaussian interstellar cloud model with the 
Doppler width, mean radial velocity and equivalent width floated 
as free parameters. A fourth free parameter allowed an additive 
shift in the continuum height if necessary. The least-squares pa- 
rameters for each sightline are shown in Table [2] and the corre- 
sponding DIB model spectrum and CH2CN velocity distribu- 
tion are shown in Figures [3] and |4j respectively. The quality of 
the fits was monitored by comparing the normalised continuum 
RMS of the observed spectra (<x c ) with the normalised RMS of 
the model fit residuals (cry). 



5.3. Hypothetical interstellar CH 2 CN distribution 

As shown by Figure[3] the model fits match the spectra to a high 
degree of accuracy for each of the eight sightlines. Variations in 
the width and structure of the DIB profile can be reproduced in 
each case, highlighted by the progression of the /18037 FWHM 
from narrow to wide in the sequence HD 186745 — » HD 183143 
— » Cyg OB2 12 — > HD 169034. The progression corresponds 
to a sequential increase in the Doppler b of the cloud model: 
16.2 -> 18.0 -> 24.0 -> 33.0 km s -1 . Such large b values would 
require unphysically large kinetic and/or turbulent gas velocities 
and clearly cannot arise in a single interstellar cloud, and this 
places some doubt on the validity of the single-cloud model used 
here. However, given the S/N of the spectra, there is no statistical 
justification for employing more than one cloud in the model. 
The Gaussian cloud model fits the DIB well as shown by cr c a 
<Tf, implying that if CH2CN is the /18037 carrier, its interstellar 
velocity distribution must be approximately Gaussian. 

Comparing the hypothetical CH2CN velocity struc- 
tures with the profiles of interstellar lines in Figure |4] 
shows that the calculated CH2CN distributions are rea- 
sonably close to the velocities of the diffuse interstellar 
clouds. Due to the lack of spectroscopic data for other 
interstellar species for five of the sightlines examined, 
Hi data were obtained from the L eiden/Argentine/Bonn 
(LAB) survey dKalberla et all I2005L data available at 



|http://www.astro.uni-bonn.de/~w ebrai/english/toolsiabsurvey.php) 
Although background H i emission contaminates these profiles, 
they at least indicate the likely velocity structure of a significant 
proportion of the neutral gas in these sightlines, aided by the 
fact that the stars are relatively distant (see Table [1). The peak 
of the calculated CH2CN distribution generally lies within c. 
10 - 15 km s of the peak Hi column density. It should be 
noted that the mean LSR velocities of the foreground clouds 
traced by Ki are positive for each sightline, which can be 
explained as a result of Galactic rotation and observational bias 
in the Galactic longitude of the targets. In Figure [4] there is 
an additional systematic small shift in the inferred CH2CN 
distribution to a higher LSR velocity as compared to the Ki. It 
is possible that this originates from the difficulty in defining the 
continuum of /18037 as it sits atop a (variable) 8040 A broader 
feature. 

Analysis of the velocity structure of interstellar gas in the ob- 
served sightlines has been performed to determine whether the 
DIB profiles can be reproduced by convolving the transitions of 
CH2CN with the observed interstellar gas distributions. The re- 
sults of the least-squares fitting (Table |2| show that a Doppler 
b parameter of ~ 20 km s _1 is required to achieve a good fit 
with the DIB. Typically, clouds traced by Ki are much narrower 
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than this ( IWeltv & Hobbsll2001l) . and even in heavily-reddened 
lines of sight such as those observed, the presence of multiple K i 
clouds separated in velocity space is insufficient to account for 
such broadening. Examining the data in left-hand panel of Figure 
|H the HD 183143 least-squares model CH 2 CN distribution 
shows gross similarities with the Tin and Can profiles. It is ev- 
ident from these spectra that Ti n has a broader and less-peaked 
distribution than the other species. In order of decreasing profile 
width and increasing structure, Ti n comes first followed by Ca n, 
then CH + , and finally Ki. This sequence is believed to reflect 
the typical conditions in the clo uds in which these species are 
generally found (see for examplelJenkinslll989[; ICrinklaw et all 
1 1 994t IWeltvl 1 1 998t IWeltv et all 120031: iPanet aul2005l) . Tin, due 
to its high condensation temperature tends to be found in the gas 
phase in hotter, more heavily shocked higher velocity clouds, 
and is depleted out onto grains in cool, quiescent diffuse clouds 
and molecular clouds. Ca has a similar condensation temperature 
to Ti and therefore comes out of the solid state to produce Ca n 
in similar energetic/shocked regions to where Tin is observed, 
but will be ionised to Ca m in more diffuse, strongly irradiated 
clouds, so it tends to be found in greatest abundances in warm 
yet moderately well-shielded regions. CH + is photodissociated 
in less dense regions, but ha s been shown to have a broader 
velocity distribution than Ki dPan et al 1 12005b . consistent with 
the formation of this molecule in war mer haloes surrounding 
cold cloud cores (ICrawfordet aUll994l) . K i is a good tracer for 
cool, quiescent, well-shielded clouds with low thermal broad- 
ening where the radiation field is attenuated and the photoioni- 
sation rate is relatively low. Based on this information and the 
breadth of the Ti n distributions shown in Figure |4] gas traced 
by Ti ii may plausibly have a velocity distribution sufficient to 
provide the Doppler broadening required for CH2CN to be the 
carrier of /18037. Tim and other gases present in regions where 
H ii is found may exist over an even greater range of velocities. 

Figure [5] shows the convolution of the 2.74 K CH 2 CN 
L Bi - X 1 A', K' a ' = transitions with the (HIRES) instrumen- 
tal PSF and with the interstellar Na i and Ti n profiles for each of 
the three sightlines. Noiseless Na i and Ti n velocity profiles were 
derived by modelling the Na i UV doublet and the Ti n 3242 and 
3384 A lines using the vapid routine dHowarth et al.L 12002). The 
accuracy of the derived Doppler b parameters was improved by 
the simultaneous modelling of two transitions of each species 
(originating in the same ground state), with differing oscillator 
strengths. Due to its similar chemistry, Nai traces the same type 
of interstellar material as K i, and was used in preference to the 
K i 7699 A lines shown in Figure [4] due to the reduced satura- 
tion and telluric contamination of the Na i UV doubl et. Atomic 
transition parameters were taken from lMortonl (120031) . 

Within the spectral signal-to-noise, the profiles in Figure [5] 
calculated by convolution of the CH2CN transitions with the 
interstellar Ti n profiles compare favourably with the diffuse in- 
terstellar band wavelength and profile, especially for HD 168 1 12 
and HD 183143. Convolution by Tin gives a significantly better 
fit than convolution by Nai in all three cases, where the nar- 
rowness of the Na i Doppler spread produces structure in the re- 
sulting contour which is not seen in the observed DIB. These 
results demonstrate that the origin band CH2CN 'Bi - X 1 A', 
K a — 1 < — transitions at 2.74 K are capable of reproduc- 
ing the narrow ,18037 DIB towards HD 183143 and HD 1681 12 
provided that the postulated CH2CN carrier molecule is dis- 
tributed in velocity space in approximately the same way as in- 
terstellar Ti ii. Due to the lack of substructure in the observed 
/18037 profile, it is very unlikely that GfbCN could instead co- 
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exist with Na i and K i in the cold neutral medium. For Cyg OB2 
8a the peak absorption wavelength and profile match is reason- 
ably good for the case of both Na i and Ti n, but there is clearly 
a significant level of structure in the calculated profiles that is 
not present in the observation. For the Ti n-convolved calcula- 
tions, the relatively small differences between the calculated and 
observed profiles could plausibly be the result of Poisson noise 
of the spectrum or other features contaminating the DIB profile 
such as other interstellar features, telluric division or flat-fielding 
residuals. 

6. Discussion 

Under the assumption that the CH2CN distribution is the same 
as that of Ti ii - with the implication that CH2CN co-exists with 
gas-phase Ti n - the fit to the observed /18037 spectrum towards 
HD 183143 and HD 1681 12 is very good (and to a lesser degree 
for Cyg OB2 8a). This result leads to two possibilities for the 
type of gas with which the postulated CH2CN might be most 
closely associated: (1) The warm, shocked 'intercloud' medium; 
for the ~ 50 interstellar clouds expect ed to be present in the 
heavily-reddened sightlines studied (see lWeltv et all 120031) . the 
warm, low density, shocked clo uds with low d e pletio ns should 
dominate the Ti n profiles {e.g. ICrinklaw et all 119941) . (2) The 
neutral ISM; Ti n is a good trac er of neutral hy drogen due to its 
ionisation potential of 13.6 eV dStokesiri978h . Without knowl- 
edge of the precise H i velocity profiles, the possibility that the 
postulated CH2CN co-exists with neutral atomic hydrogen gas 
cannot be ruled out. 

There are no reports of past attempts to detect CH2CN 
or CH 2 CN in the diffuse ISM. The neutral radical CH 2 CN 
has been o bserved in the dense TMC-1 and Sgr B 2 molecular 
clouds (see llrvine et all 119881: iTurner et all 1 19901) through its 
pure rotational transitions. Fractional abundances are ~ 10~ 9 to 
10" 11 relative to Hz, and CH2CN was found to be similarly dis- 
tributed to CH 3 CN and C 4 H. The synthesis of CH 2 CN in the 
shocked diffuse interstellar clouds traced by Ti 11 could be driven 
by the release of hydrocarbons into the gas phase during su- 
pernova shock-induced collisions between carbon aceous grains 
dDulev & Williamsl! 1984 [Hall & Wimamslll995l) . Given a suf- 
ficient abundance of C2H4 in the gas phase, the reaction 

C 2 Hj; + N — > CH 3 CN + + H, (1) 

with a rate coefficient of ~ 10 -10 cm 3 s _1 dHerbst & Leungl 

119901) . could proceed rapidly enough to produce sufficient quan- 
tities of CH3CN + to allow dissociative recombination (DR) 

CH 3 CN + + e — CH 2 CN + H (2) 

to occur. This mechanism for CH 2 CN production is only vi- 
able if single H-atom loss is a significant channel in the DR of 
CH 3 CN + (a process which is yet to be studied in detail), and 
if recombination occurs before CH 3 CN + is destroyed by react- 
ing with neutral hydrogen. The dominant mechanisms by which 
CH2CN and CH 3 CN are expected to be destroyed in the diffuse 
ISM are by photodissociation and by reaction with C + . 

The most obvious pathway for production of the 
cyanomethyl anion CH2CN is by radiative electron attachment 
to CH 2 CN: 

CH 2 CN + e — > CH 2 CN + hv, (3) 

though an alternative route exists via dissociative attachment 
to CH 3 CN: 

CH 3 CN + e — -» CH 2 CN + H. (4) 
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The last step was studied bv lSailer et al.l d2003l) in the labora- 
tory. The dissociative electron attachment cross-section peaked 
at 4 x 10~ 23 m 2 for electron energies of 3.2 eV, which provides 
an efficient route to the formation of CH2CN from CH3CN. 
The cyanomethyl anion production channel is ~ 500 times 
more efficient than those of the other possible anionic fragments 
(CHCN , CCN , CN and CH 3 ) for electrons energies between 
1 and 10 eV. In the interstellar context, photon and cosmic ray 
impact on grains could provide a source of electrons with a few 
eV energy. 

If the dominant chemical reactions involving CH2CN are 
assumed to be formation of the anion via radiative attachment 
(Equation (O) and destruction via photodetachment, then in 
equilibrium, the anion-to-neutral ratio is given by 

[CH 2 CN-] = an, 

[CH 2 CN] r 

where a is the rate of radiative electron attachment, n e is the 
electron density and F is the photodetachment rate. Following 
s-wave el ectron capture by CH2C N into an excited state of 
CH2CN dTerzieva & Herbstl l2000h . if subsequent radiative re- 
laxation occurs rapidly, either by vibrational transitions of the vi- 
brationally excited ground state or by electronic transitions from 
the dipole-bound state to the ground state, then a can be esti- 
mated using Equation (6) of iTerzieva & Herbstl d2000h . In the 
diffuse ISM (at T = 100 K), this yields a value of a = 2.2 x 10~ 7 
cirT 3 s , assuming that the exited anion undergoes radiative 
stabilisation before the electron detaches back to the contin- 
uum. The photodetachment rate F has been calculated to be 
approximately 1.2 x 10~ 8 s -1 (E. Herbst & T. J. Millar, pri- 
vate communi cation) in the standard interstellar radiation field 
dDraineL [1978). Both parameters are uncertain due to the lack 
of laboratory or theoretical data for this molecule, but neverthe- 
less permit a rough estimate of the fraction of CH2CN likely 
to occur in anionic form. Using an electro n density in t he dif - 
fuse ISM of between 0.04 and 0.23 cm 3 dWeltv et al.L 120031) . 
Equation (O yields an anion-to-neutral ratio between 0.7 and 
4.2. If CH2CN is present in the warm n eutral medium trac ed by 
Tin (and Can) where T ~ 4000 K (see lWeltv et"ailll996l) then 
a = 3.4x 10 -8 cm -3 s _1 and the calculated anion-to-neutral ratio 
is 0.1-0.7. 

A further test of the plausibility of CH2CN as carrier of 
the /18037 DIB is whether a feasible CH2CN abundance can 
be deduced from the equivalent width of the observed DIB. 
Polyatomic molecules are likely to be present in all of the 
heavily-reddened sightlines towards which /i8037 was observed. 
Cy g OB 12, Cyg OB 2 5 and HP 183143 have been studied 
by iMcCall et all d!998l l2002al) . particularly with reference to 
interstellar H3 + , detected in all three sightlines with Doppler 
b ~ 10 km s _1 which is substantially broader than the 13 CO and 
HCO + lines observed towards Cyg O B2 12. The HCO + col umn 
density was found to be ~ 10" cirr 2 . IScappini et all d2002l) de- 
termined a likely total number density of n > 10 4 cm -3 for the 
gas containing the CO and HCO + , with the inference of dense, 
compact molecular clumps within a more diffuse medium. The 
approximate molecular hydrogen c olumn density towards Cyg 
OB2 12 (from lMcCall"etail l2002al) is 6.5 x 10 21 cirT 2 such that 
the fractional HCO + abundance is ~ 1.5 x 10 _11 nH 2 - 

The equivalent width of the ,18037 DIB is 38.3 mA (see 
Table which requires a column density of Af(CH 2 CN~) = 
6.7 x 10 1() cm -2 // where / is the oscillator strength of the tran- 
sition. Transitions between the gr ound and dipole-bo und states 
are strongly observed in the lab dLvkke et all 1 1 987b . suggest- 
ing that the value of / should be relatively large. Therefore, if 



/ = 0.5 and assuming a CH2CN anion-to-neutral ratio of 1, a 
molecular fraction of CH2CN of 2 x 10~ h «h 2 1S required to- 
wards Cyg OB2 12 to produce the observed /18037 equivalent 
width. Thus, provided the 1 B 1 - X 'A' transition is strong, the 
required interstellar CH2CN fractional column density is com- 
parable to that of HCO + observed in the Cyg OB2 12 sightline. 
Averaged along the sightline, the required interstellar CH2CN 
fractional abundance is at the lower end of the range of val- 
ues found in dense molecular clouds (~ 10~ n - 10~ 9 hh 2 )- If 
diffuse and dense cloud molecular abundances are similar, then 
enough CH2CN can probably be produced to create the /18037 
DIB along heavily-reddened sightlines. 

The spectra observed bv lMcCall et al.l d2002al) show clearly 
that the interstellar material traced by H3 has a broader veloc- 
ity distribution than the molecules (CO, C2, CN) which typi- 
cally trace dense molecular material. The Hi line FWHM ob- 
served towards Cyg OB2 12, Cyg OB2 5 and HD 183143 range 
from 8 to 15 km s" 1 whereas the CN line FWHM are all less 
than 4 km s~'. lMcCall et al.l d2002al) conclude that 'the chem- 
istry that leads to H3 is completely decoupled from that which 
is responsible for these heavier diatomics'. These factors lend 
support to the possibility that other molecules may be present in 
significant abundances in the more diffuse regions surrounding 
conventional molecular clouds. 

Very rapid production mechanisms are required for 
molecules such as CH2CN and its anion to exist in the warm 
diffuse regions traced by Ti n, though PAHs and related large or- 
ganic species will survive for much longer in the UV radiation 
field with the potential to give rise to diffuse interstellar bands. 
The interst ellar media containing Ti n and Ca n are largely co- 
spatial (see lCrinklaw et al.L [19941 J Albert et al.L 1 1993b . and these 
refractor y ions are observe d predominantly in the warm neu- 
tral ISM. IWeltv et al.l d 19961) found the mean upper temperature 
limit of the clouds in their Can survey to be 4100 K, though 
temperatures as low as a few hundred K were found in many 
cases. Titanium and calcium arise in the g as phase predomi- 
nantly due to grain destruction dStokesL[l978l) . and it is plausible 
that molecules and DIB carriers may also arise in this way. 

There is no evidence in the literature for correl ations between 
DIB strengths and Ti 11 column densities, thoug h lHerbiddl993h 
found that the 5780 and 5797 DIBs do not correlate with the 
level of titanium depletion. This study suggests that further in- 
vestigation of a possible link between DIBs, molecules and Tin 
is warranted. 



7. Conclusion 

Using the rotational constants of iLvkke et al.1 dl987l) . the ab- 
sorption spectrum arising from the 'Bi - X 'A' transition of 
CH2CN has been computed for CH2CN in equilibrium with 
the CMB at 2.74 K. The peak of the intrinsic absorption profile 
(at 8037.78 A) is consistent with the peak of the diffuse interstel- 
lar absorption feature at 8037.8 + 0.15 A found in the spectra of 
eight heavily-reddened stars. If CH2CN occurs in the ISM with 
the statistical ortho : para abundance ratio of 3 : lthe'Bi-X'A' 
transitions with K'J = 1 would produce strong spectral features 
at wavelengths of approximately 8024.8 and 8049.6 A that are 
not present in the observed interstellar spectra. However, if the 
K'J = 1 level is populated based on a Boltzmann distribution 
at 2.74 K, then the K a = < — 1 and K a = 2 < — 1 sub- 
bands are sufficiently weak to be compatible with their absence 
from the observed spectra. A thermal (Boltzmann) distribution 
of K'J levels would be expected if an efficient mechanism exists 
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for conversion between the ortho and para forms of CH2CN . 
Alternatively, the formation mechanism(s) for CTTCN may in- 
troduce a skewed (i.e. non 3:1) ratio of ortho and para forms. 
Higher signal-to-noise (> 2000) spectroscopic observations of 
the /18037 region will assist in the search for the K'J = 1 sub- 
bands. Without detailed laboratory or theoretical studies, it is 
unclear whether or not the CH>CN formation mechanism(s) 
would result in the required deviation of the ortho-to-para ratio 
from the statistical value of 3 : 1 . 

The cyanomethyl anion CH2CN is a plausible candidate 
for the carrier of the narrow diffuse interstellar band located at 
8037 .8 A, provided a spectral line-broadening mechanism can be 
accounted for that produces an approximately Gaussian broad- 
ening corresponding to a Dopplerfr of between 16 and 33 km s _1 
which corresponds with the observed Ti 11 and Ca 11 profiles to- 
wards HD 183143 and HD 1681 12, and is comparable to the H+ 
linewidths towards Cyg OB2 12, Cyg OB2 5 and HD 183143. 

Measurement of the oscillator strength of the 'Bi - X 'A' 
transition would enable the calculation of the interstellar col- 
umn density of CH2CN that would be required to produce the 
/18037 DIB. Assuming the transition is strong (/ ~ 0.5), and 
the CH2CN anion-to-neutral ratio is large (~ 1, as suggested by 
preliminary calculations of the anion formation and destruction 
rates), then towards Cyg OB2 12 the required average fractional 
abundance of neutral CH2CN relative to Hi is 2 x 10 , which 
falls at the lower end of the observed range in dense clouds. 

The spectroscopic consistency between the ,18037 DIB and 
the calculated transitions of CH2CN cannot be considered an 
assignment, but provides strong motivation for further study of 
this molecular anion including a search for the pure rotational 
transitions of interstellar CH2CN . 
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